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The supersymmetric extension to the Standard Model offers a promising cold dark 
matter candidate, the lightest neutralino. f will review the prospects for the de- 
tection of this candidate in both accelerator and direct detection searches. 



1. Introduction 

Although there are many reasons for considering supersymmetry as a can- 
didate extension to the standard model of strong, weak and electromagnetic 
interactions 1 , one of the most compelling is its role in understanding the 
hierarchy problem 2 namely, why/how is my/ <C Mp. One might think 
naively that it would be sufficient to set my/ <C Mp by hand. However, 
radiative corrections tend to destroy this hierarchy. For example, one-loop 
diagrams generate 

Sm 2 w = O A 2 » my/ (1) 

where A is a cut-off representing the appearance of new physics, and the 
inequality in (1) applies if A ~ 10 3 TeV, and even more so if A ~ mauT ~ 
10 16 GeV or ~ Mp ~ 10 19 GeV. If the radiative corrections to a physical 
quantity are much larger than its measured values, obtaining the latter 
requires strong cancellations, which in general require fine tuning of the 
bare input parameters. However, the necessary cancellations are natural in 
supersymmetry, where one has equal numbers of bosons B and fermions F 
with equal couplings, so that (1) is replaced by 

8m 2 w = 0(^) \ml-m 2 F \ . (2) 
The residual radiative correction is naturally small if \m 2 B — m 2 F \ < 1 TeV 2 . 
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In order to justify the absence of interactions which can be responsible 
for extremely rapid proton decay, it is common in the minimal supersym- 
metric standard model (MSSM) to assume the conservation of R-parity. If 
R-parity, which distinguishes between "normal" matter and the supersym- 
metric partners and can be defined in terms of baryon, lepton and spin as 
R = (_i) 3B + i + 2 ' s j i s unbroken, there is at least one supersymmetric par- 
ticle (the lightest supersymmetric particle or LSP) which must be stable. 
Thus, the minimal model contains the fewest number of new particles and 
interactions necessary to make a consistent theory. 

There are very strong constraints, however, forbidding the existence of 
stable or long lived particles which are not color and electrically neutral 3 . 
Strong and electromagnctically interacting LSPs would become bound with 
normal matter forming anomalously heavy isotopes. Indeed, there are very 
strong upper limits on the abundances, relative to hydrogen, of nuclear 
isotopes 4 , n/n H < 1CT 15 to 1(T 29 for 1 GcV < m < 1 TcV. A strongly 
interacting stable relic is expected to have an abundance n/nn Si 10~ 10 
with a higher abundance for charged particles. 

There are relatively few supersymmetric candidates which are not col- 
ored and are electrically neutral. The sneutrino 5 is one possibility, but in 
the MSSM, it has been excluded as a dark matter candidate by direct 6 and 
indirect 7 searches. In fact, one can set an accelerator based limit on the 
sneutrino mass from neutrino counting, > 44.7 GeV 8 . In this case, the 
direct relic searches in underground low-background experiments require 
mo > 20 TeV 6 . Another possibility is the gravitino which is probably the 
most difficult to exclude. I will concentrate on the remaining possibility in 
the MSSM, namely the neutralinos. 

2. Parameters 

The most general version of the MSSM, despite its minimality in particles 
and interactions contains well over a hundred new parameters. The study 
of such a model would be untenable were it not for some (well motivated) 
assumptions. These have to do with the parameters associated with super- 
symmetry breaking. It is often assumed that, at some unification scale, all 
of the gaugino masses receive a common mass, m^. The gaugino masses 
at the weak scale are determined by running a set of renormalization group 
equations. Similarly, one often assumes that all scalars receive a common 
mass, mo, at the GUT scale. These too are run down to the weak scale. 
The remaining parameters of importance involve the Higgs sector. There is 
the Higgs mixing mass parameter, fi, and since there are two Higgs doublets 
in the MSSM, there are two vacuum expectation values. One combination 
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of these is related to the Z mass, and therefore is not a free parameter, 
while the other combination, the ratio of the two vevs, tan/3, is free. 

If the supersymmetry breaking Higgs soft masses are also unified at the 
GUT scale (and take the common value mo), then /j, and the physical Higgs 
masses at the weak scale are determined by electroweak vacuum conditions. 
This scenario is often referred to as the constrained MSSM or CMSSM. 
Once these parameters are set, the entire spectrum of sparticle masses at 
the weak scale can be calculated. 



3. Neutralinos 



There are four neutralinos, each of which is a linear combination of the 
R = — 1 neutral fermions, 3 : the wino W 3 , the partner of the 3rd component 
of the SU(2)l gauge boson; the bino, B, the partner of the U(l)y gauge 
boson; and the two neutral Higgsinos, Hi and H 2 . Assuming gaugino 
mass universality at the GUT scale, the identity and mass of the LSP are 
determined by the gaugino mass m 1 / 2 , ^, and tan/3. In general, neutralinos 
can be expressed as a linear combination 



X = aB + (3W 3 + 7 ifi + 6H 2 



(3) 



The solution for the coefficients a, /3, 7 and (5 for neutralinos that make up 
the LSP can be found by diagonalizing the mass matrix 
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where Mi(M2) is a soft supersymmetry breaking term giving mass to the 
U(l) (SU(2)) gaugino(s). In a unified theory Mi = M 2 at the unification 
scale (at the weak scale, Mi ~ |^-M 2 ). As one can see, the coefficients 
a,/3,7, and S depend only on m\/ 2 , /U, and tan/3. 

In Figure 1 9 , regions in the M 2 ,fi plane with tan/3 = 2 are shown 
in which the LSP is one of several nearly pure states, the photino, 7, the 
bino, B, a symmetric combination of the Higgsinos, ff(i2), ° r the Higgsino, 
S = sin (3H\ + cos (3H 2 ■ The dashed lines show the LSP mass contours. The 
cross hatched regions correspond to parameters giving a chargino {W ± 1 H^) 
state with mass < 45GeU and as such are excluded by LEP 10 . This 
constraint has been extended by LEP 11 and is shown by the light shaded 
region and corresponds to regions where the chargino mass is < 104 GeV. 
The newer limit does not extend deep into the Higgsino region because 
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of the degeneracy between the chargino and neutralino. Notice that the 
parameter space is dominated by the B or Hi 2 pure states and that the 
photino only occupies a small fraction of the parameter space, as docs 
the Higgsino combination S. Both of these light states are experimentally 
excluded. 




Figure 1. Mass contours and composition of nearly pure LSP states in the MSSM 9 . 



4. The Relic Density 

The relic abundance of LSP's is determined by solving the Boltzmann 
equation for the LSP number density in an expanding Universe. The 
technique 12 used is similar to that for computing the relic abundance of 
massive neutrinos 13 . The relic density depends on additional parameters in 
the MSSM beyond mi/ 2 ,/U, an d tan/3. These include the sfermion masses, 
m j and the Higgs pseudo-scalar mass, m,A a , derived from m (and my 2 ). 

a In general, the relic density depends on the supcrsymmctry-breaking tri-linear masses 
A (also assumed to be unified at the GUT scale) as well as two phases M and 8 a- 
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To determine the relic density it is necessary to obtain the general anni- 
hilation cross-section for neutralinos. In much of the parameter space of 
interest, the LSP is a bino and the annihilation proceeds mainly through 
sfermion exchange. Because of the p-wave suppression associated with Ma- 
jorana fermions, the s-wave part of the annihilation cross-section is sup- 
pressed by the outgoing fermion masses. This means that it is necessary 
to expand the cross-section to include p-wave corrections which can be 
expressed as a term proportional to the temperature if neutralinos are in 
equilibrium. Unless the neutralino mass happens to lie near near a pole, 
such as m x ~ mz/2 or mh/2, in which case there are large contributions 
to the annihilation through direct s-channel resonance exchange, the domi- 
nant contribution to the BB annihilation cross section comes from crossed 
i-channel sfermion exchange. 

Annihilations in the early Universe continue until the annihilation rate 
T ~ <rvn x drops below the expansion rate given by the Hubble param- 
eter, H. For particles which annihilate through approximate weak scale 
interactions, this occurs when T ~ m x /20. Subsequently, the relic density 
of neutralinos is fixed relative to the number of relativistic particles. As 
noted above, the number density of neutralinos is tracked by a Boltzmann- 
like equation, 

Tt =- 3 |"-M(^ 2 -^o) (5) 

where no is the equilibrium number density of neutralinos. By defining the 
quantity / = n/T 3 , we can rewrite this equation in terms of x, as 

I -™»(^) " 2 (f'-fS) (0) 

The solution to this equation at late times (small x) yields a constant value 
of /, so that n oc T 3 . The final relic density expressed as a fraction of the 
critical energy density can be written as 3 

n x h 2 ~ 1.9 x io- 11 (^X n)' 2 ( Ge ^ \ (7) 

x \Tj 1 \ax f + \bx) ) 1 ; 

where (T x /T 7 ) 3 accounts for the subsequent reheating of the photon tem- 
perature with respect to \, due to the annihilations of particles with mass 
m < Xfm x 14 . The subscript / refers to values at freeze-out, i.e., when 
annihilations cease. The coefficients a and b are related to the partial wave 
expansion of the cross-section, <jv = a + bx + . . .. Eq. (7 ) results in a very 
good approximation to the relic density expect near s-channel annihilation 
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poles, thresholds and in regions where the LSP is nearly degenerate with 
the next lightest supersymmetric particle 15 . 

5. Phenomenological and Cosmological Constraints 

For the cosmological limits on the relic density I will assume 

0.1 < n x h 2 < 0.3. (8) 

The upper limit being a conservative bound based only on the lower limit 
to the age of the Universe of 12 Gyr. Indeed, most analyses indicate that 
^matter Ss 0.4 — 0.5 and thus it is very likely that fl x h 2 < 0.2. One should 
note that smaller values of il x h 2 are allowed, since it is quite possible that 
some of the cold dark matter might not consist of LSPs. 

The calculated relic density is found to have a relevant cosmological 
density over a wide range of susy parameters. For all values of tan/3, 
there is a 'bulk' region with relatively low values of mi/2 and mo where 
0.1 < Q x h 2 < 0.3. However there are a number of regions at large values 
of mi/2 and/or m where the relic density is still compatible with the 
cosmological constraints. At large values of m 1 / 2 , the lighter stau, becomes 
nearly degenerate with the neutralino and co-annihilations between these 
particles must be taken into account 16 . For non-zero values of Aq, there 
are new regions for which \ ~ t coannihilations are important 17 . At large 
tan /?, as one increases mi/ 2 , the pseudo-scalar mass, tha begins to drop 
so that there is a wide funnel- like region (at all values of m ) such that 
2m x « uia and s-channel annihilations become important 18,19 . Finally, 
there is a region at very high mo where the value of \i begins to fall and 
the LSP becomes more Higgsino-like. This is known as the 'focus point' 
region 20 . 

As an aid to the assessment of the prospects for detecting sparticles at 
different accelerators, benchmark sets of supersymmetric parameters have 
often been found useful, since they provide a focus for concentrated discus- 
sion. A set of proposed post-LEP benchmark scenarios 21 in the CMSSM 
are illustrated schematically in Fig. 2. Five of the chosen points are in the 
'bulk' region at small m\/2 and mo, four are spread along the coannihila- 
tion 'tail' at larger mi/2 for various values of tan/3. This tail runs along 
the shaded region in the lower right corner where the stau is the LSP and 
is therefore excluded by the constraints against charged dark matter. Two 
points arc in rapid-annihilation 'funnels' at large TO1/2 and mo- At large 
values of mo, the focus-point region runs along the boundary where elec- 
troweak symmetry no longer occurs (shown in Fig. 2 as the shaded region 
in the upper left corner). Two points were chosen in the focus-point region 
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Figure 2. Schematic overview of the CMSSM benchmark points proposed in . The 
points are intended to illustrate the range of available possibilities. The labels correspond 
to the approximate positions of the benchmark points in the (m 1 ^2t m o) plane. They 
also span values of tan j3 from 5 to 50 and include points with fi < 0. 

at large too- The proposed points range over the allowed values of tan/3 
between 5 and 50. The light shaded region corresponds to the portion of 
parameter space where the relic density Q x h 2 is between 0.1 and 0.3. 

The most important phenomenological constraints are also shown 
schematically in Figure 2. These include the constraint provided by the 
LEP lower limit on the Higgs mass: ma > 114.1 GeV 22 . This holds in 
the Standard Model, for the lightest Higgs boson h in the general MSSM 
for tan/3 < 8, and almost always in the CMSSM for all tan/3. Since mh 
is sensitive to sparticle masses, particularly mj, via loop corrections, the 
Higgs limit also imposes important constraints on the CMSSM parameters, 
principally m\i2 as seen by the dashed curve in Fig. 2. The constraint 
imposed by measurements of b — > 57 23 also exclude small values of m\ii- 
These measurements agree with the Standard Model, and therefore pro- 
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Figure 3. Compilation of phenomenological constraints on t he CMSSM for tan j3 



10, fi > 0, assuming Aq = 0,mt = 175 GeV and mj(mj) 



MS 



4.25 GeV. The near- 



vertical lines are the LEP limits m„± = 103.5 GeV (clashed and black) 11 , and = 
114.1 GeV (dotted and red) 22 . Also, in the lower left corner we show the mg = 99 GeV 
contour 28 . In the dark (brick red) shaded regions, the LSP is the charged fi, so this 
region is excluded. The light (turquoise) shaded areas are the cosmologically preferred 
regions with 0.1 < f!h 2 < 0.3 19 . The medium (dark green) shaded regions are excluded 
by b — > S7 23 . The shaded (pink) region in the upper right delineates the 2 a range of 
9m -2- 



vide bounds on MSSM particles, such as the chargino and charged Higgs 
masses, in particular. Typically, the b — > sj constraint is more impor- 
tant for fi < 0, but it is also relevant for /i > 0, particularly when tan/3 
is large. The BNL E821 experiment reported last year a new measure- 
ment of = \{g^ — 2) which deviated by 2.6 standard deviations from 
the best Standard Model prediction available at that time 24 . The largest 
contribution to the errors in the comparison with theory was thought to 
be the statistical error of the experiment, which has been significantly re- 
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duced just recently 25 . However, it has recently been realized that the sign 
of the most important pseudoscalar-meson pole part of the light- by- light 
scattering contribution 26 to the Standard Model prediction should be re- 
versed, which reduces the apparent experimental discrepancy to about 1.6 
standard deviations (Sa^ x 10 10 = 26 ± 16). With the new data, the dis- 
crepancy with theory ranges from 1.6 to 2.6 a, i.e., Sa^ x 10 10 = 26 ± 10 
to 17 ± ll 25 . This constraint excludes very small values of and ma- 
in Fig. 2, the g — 2 constraint is shown schematically by the dotted line. 
It may also exclude very large values of the parameters as well as negative 
values of [i, if the discrepancy holds up. 

Following a previous analysis 27 , in Figure 3 the TO1/2 — mo parame- 
ter space is shown for tan/3 = 10. The dark shaded region (in the lower 
right) corresponds to the parameters where the LSP is not a neutralino but 
rather a f^. The cosmologically interesting region at the left of the figure 
is due to the appearance of pole effects. There, the LSP can annihilate 
through s-channel Z and h (the light Higgs) exchange, thereby allowing a 
very large value of m . However, this region is excluded by phenomeno- 
logical constraints. Here one can see clearly the coannihilation tail which 
extends towards large values of m 1 / 2 - In addition to the phenomenological 
constraints discussed above, Figure 3 also shows the current experimental 
constraints on the CMSSM parameter space due to the limit m x ± > 103.5 
GeV provided by chargino searches at LEP 11 . LEP has also provided lower 
limits on slepton masses, of which the strongest is mg > 99 GeV 28 . This 
is shown by dot-dashed curve in the lower left corner of Fig. 3. 

As one can see, one of the most important phenomenological constraint 
at this value of tan (i is due to the Higgs mass (shown by the nearly verti- 
cal dot-dashed curve). The theoretical Higgs masses were evaluated using 
FeynHiggs 29 , which is estimated to have a residual uncertainty of a couple 
of GeV in to/,. The region excluded by the b — > sj constraint is the dark 
shaded (green) region to the left of the plot. 

As many authors have pointed out 30 , a discrepancy between theory and 
the BNL experiment could well be explained by supersymmetry. As seen 
in Fig. 3, this is particularly easy if > 0. The medium (pink) shaded 
region in the figure corresponds to the overall allowed region by the new 
experimental result: —5 < 5a^ x 10 10 < 46. The two solid lines within the 
shaded region corresponds to the central values <5a M x 10 10 = 17 and 26 
respectively. The optimistic 2a lower bound of <5a M x 10 10 = 6 is shown as 
a dashed curve. 

As discussed above, another mechanism for extending the allowed 
CMSSM region to large m x is rapid annihilation via a direct-channel pole 
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Figure 4. As in Fig. 3 for tan f} = 50. 



when m x ~ ^m^ 18 > 19 . This may yield a 'funnel' extending to large m-i/2 
and mo at large tan /?, as seen in Fig. 4. 

In principle the true input parameters in the CMSSM are: fj,,mi,rri2, 
and B, where mi and m,2 are the Higgs soft masses (in the CMSSM 
mi = m.2 = mo and B is the susy breaking bilinear mass term). In this 
case, the electroweak symmetry breaking conditions lead to a prediction of 
Mz-i tan (3 ,and tua- Since we are not really interested in predicting Mz, 
it is more useful to assume instead the following CMSSM input param- 
eters: Mz,mi,m2, and tan/3 again with mi = ni2 = Wo- I n this case, 
one predicts fi,B, and tua- However, one can generalize the CMSSM case 
to include non-universal Higgs masses (NUHM), in which case the input 
parameters become:Mz, /x, tua, and tan/3 and one predicts mi, m-2, and B. 

The NUHM parameter space was recently analyzed 31 and a sample of 
the results found is shown in Fig. 5. While much of the cosmologically 



11 



preferred area with /i < is excluded, there is a significant enhancement in 
the allowed parameter space for /i > 0. 



tan p = 10 , m 1/2 = 300 , m = 100 




jj, (GeV) 

Figure 5. Compilations of phenomenological constraints on the MSSM with NUHM 
in the (fi, vtia) plane for tan j3 = 10 and mo = 100 GeV, mi/j = 300 GeV, assuming 

Ao = 0, mt = 175 GeV and fii,(?7lj,)gj{| = 4.25 GeV. The shading is as described in 
Fig. 3. The (blue) solid line is the contour m x = near which rapid direct-channel 

annihilation suppresses the relic density. The dark (black) dot-dashed line indicates when 
one or another Higgs mass-squared becomes negative at the GUT scale: only lower \fi\ 
and larger rriA values are allowed. The crosses denote the values of fi and found in 
the CMSSM. 



5.1. Detection 

Because the LSP as dark matter is present locally, there are many avenues 
for pursuing dark matter detection. Direct detection techniques rely on 
an ample neutralino-nucleon scattering cross-section. The effective four- 
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fermion lagrangian can be written as 

+ aziXXmi + ol^xi^X^H 

+ "5 l XX%7 5 9* + a6iXl 5 XQiqi (9) 

However, the terms involving an, an, a$i, and a§i lead to velocity depen- 
dent elastic cross sections. The remaining terms are: the spin dependent 
coefficient, a 2 i and the scalar coefficient a^i- Contributions to a 2 i are 
predominantly through light squark exchange. This is the dominant chan- 
nel for binos. Scattering also occurs through Z exchange but this channel 
requires a strong Higgsino component. Contributions to a$i are also dom- 
inated by light squark exchange but Higgs exchange is non-negligible in 
most cases. 




m% [GeV] 

Figure 6. Limits from the CDMS 34 and Edelweiss 35 experiments on the neutralino- 
proton elastic scattering cross section as a function of the neutralino mass. The Edelweiss 
limit is stronger at higher m x . These results nearly exclude the shaded region observed 
by DAMA 36 . The theoretical predictions lie at lower values of the cross section. 
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The results from a CMSSM and MSSM analysis 32 ' 33 for tan /? = 3 and 
10 are compared with the most recent CDMS 34 and Edelweiss 35 bounds in 
Fig. 6. These results have nearly entirely excluded the region purported 
by the DAMA 36 experiment. The CMSSM prediction 32 is shown by the 
dark shaded region, while the NUHM case 33 is shown by the larger lighter 
shaded region. 

I conclude by showing the prospects for direct detection for the 
benchmark points discussed above 37 . Fig. 7 shows rates for the elastic 
spin-independent scattering of supersymmetric relics, including the pro- 
jected sensitivities for CDMS II 38 and CRESST 39 (solid) and GENIUS 40 
(dashed). Also shown are the cross sections calculated in the proposed 
benchmark scenarios discussed in the previous section, which are consid- 
erably below the DAMA 36 range (10~ 5 — 10~ 6 pb). Indirect searches for 
supersymmetric dark matter via the products of annihilations in the galac- 
tic halo or inside the Sun also have prospects in some of the benchmark 
scenarios 37 . 




Figure 7. Elastic spin-independent scattering of supersymmetric relics on protons cal- 
culated in benchmark scenarios 37 , compared with the projected sensitivities for CDMS 
II 38 and CRESST 39 (solid) and GENIUS 40 (dashed). The predictions of our code (blue 
crosses) and Neutdriver 41 (red circles) for neutralino-nucleon scattering are compared. 
The labels A, B, ...,L correspond to the benchmark points as shown in Fig. 2. 
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